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Abstract. Novel techniques are presented, which identify the chiral power-counting regime (PCR), 
i and realize the existence of an intrinsic energy scale embedded in lattice QCD results that extend 

outside the PCR. The nucleon mass is considered as a benchmark for illustrating this new approach. 
Using finite-range regularization, an optimal regularization scale can be extracted from lattice 
CsJ ■ simulation results by analyzing the renormalization of the low energy coefficients. The optimal scale 

\ allows a description of lattice simulation results that extend beyond the PCR by quantifying and thus 

handling any scheme-dependence. Preliminary results for the nucleon magnetic moment are also 
examined, and a consistent optimal regularization scale is obtained. This indicates the existence of 
an intrinsic scale corresponding to the finite size of the source of the pion cloud. 
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cn INTRODUCTION 

> ' 

^ . Chiral perturbation theory (xPT) provides a formal approach to counting the powers of 

[ low energy momenta and quark masses, so that an ordered expansion in the quark mass 

^ I (or pion mass, through lUq <x m|) can be constructed. ;tPT indicates that, in general, the 

■ most singular non-analytic contributions to hadron properties lie in the low order meson 
loop of the hadron. For example, the leading non-analytic behaviour of a baryon mass is 

3/2 T 

proportional to niq or m^. More generally, baryon masses can be written as an ordered 
expansion in quark mass or m|. 
^ I To establish a model-independent framework in ;tPT, the expansion must display the 

^ ■ properties of a convergent series for the terms considered. There is a power-counting 

- - -' regime (PCR) where the quark mass is sufficiently small such that higher order terms in 

the expansion are negligible beyond the order calculated. That is, within the PCR, the 
truncation of the chiral expansion is reliable to a prescribed precision [1, 2]. 

The asymptotic nature of the chiral expansion places the focus on the first few terms 
of the expansion. In chiral effective field theory (;^EFT), calculations beyond one-loop 
are time-consuming and there are no two-loop calculations which incorporate the effects 
of placing a baryon in a finite volume. With only a few terms of the expansion known for 
certain, knowledge of the PCR of ;tEFT is as important as knowledge of the expansion 
itself. 

Numerical simulations of QCD on a spacetime lattice are complemented by ;t;EFT 
through a model-independent formalism for connecting lattice simulation results to the 
physical region. Simulations at finite volume and a variety of quark masses are related 



to the infinite- volume and physical quark masses through this formalism. However, it 
is accurate only within the PGR of the truncated expansion. In particular, truncated 
expansions are typically applied to a wide range of quark (or pion) masses with little 
regard to a careful determination of the PGR. 

Through the consideration of a variety of renormalization schemes and associated pa- 
rameters, new techniques to identify the PGR have been established in Ref [3]. Gonsider 
the nucleon mass: a benchmark for illustrating this approach. 



EFFECTIVE FIELD THEORY FOR NUCLEONS 



Using the standard Gell-Mann-Oakes-Renner relation connecting quark and pion 



masses, niq oc m| [4], the formal chiral expansion of the nucleon mass can be written as 



the sum of a polynomial expansion in m| and the meson-loop integral contributions: 

Mn = {ao + a2ml + a4m^ + ^(m^) } + + ^^ + ^tad ■ (1) 

The pion cloud corrections are considered in the heavy-baryon limit, with loop in- 
tegrals, Lat, La and 'Lfad- The coefficients a; of the analytic polynomial in Eq. (1) are 
related to the low energy constants of ;fPT. In this investigation, they will be deter- 
mined by fitting to lattice simulation results. These coefficients will be referred to as 
the residual series coefficients. These bare coefficients undergo renormalization due to 
contributions from the loop integrals Lat, La and 'Ltad- 

Under the most general considerations, each loop integral, when evaluated, produces 
an analytic polynomial in m| and non-analytic terms. Using finite-range regularization, 
expanding out each integral yields: 

(3) 

^ta, = C2ml [xt ^ / d'^^^^^lP - ^0 = C2{b'2ml + bWn + Xtm^ log + ^{m^ , 

(4) 

where u{k^;A) is a regulator function described below. These integrals are expressed in 
terms of the pion energy (o{k) — ^JW+n^ and the delta-nucleon mass splitting A. Xn, 
Xi^ and Xt denote the model-independent chiral coefficients of the terms that are non- 
analytic in the quark mass. The hi coefficients are renormalization-scheme dependent, 
as are the a, coefficients. 

The process of renormalization in finite-range regularized ;tEFT proceeds by com- 
bining the renormalization-scheme dependent coefficients to provide the physical low 
energy coefficients, which are denoted as c/. Thus, the nucleon mass expansion takes on 



the standard forai: 



= Co + C2ml + XNml + C4m^ + ^- -^X^ + ^2;^?^ log ^ + ^(m^) . (5) 

By comparing Eqs. (1) and (5), the following renormalization procedure is obtained: 

CO = UQ + b^ + b^, (6) 
C2 = a2 + b'^2+b^ + C2b'2, (7) 
C4 = aA + yi + b^ + C2b\, Qic. (8) 

The coefficients c; are scheme-independent quantities, and this property will be demon- 
strated when determined within the PCR. 

For the purposes of this investigation, the regulator function u{k;A) is chosen to be a 
dipole form: 



M(fc;A)= (^1 + ^j . (9) 

More extensive discussions of different finite-range regulator functional forms can be 
found in Refs [1, 5, 6, 7]. An analysis of results for the nucleon mass for a family of 
dipole-like forms is reported in Ref [3]. 

In order to accommodate the effect of the finite volume, the continuous loop integrals 
are transformed into a sum over discrete momentum values ki = niln/L, for integers 
Hi. The difference between a loop sum and its corresponding loop integral is the finite- 
volume correction, which should vanish for all integrals as m^^L becomes large [8]. 



THE INTRINSIC SCALE 

The xEFI extrapolation scheme will be used in conjunction with lattice QCD data from 
CP-PACS [9] and PACS-CS [10] to predict the nucleon mass Mj^ for any value of m%. 
When using lattice QCD results beyond the PCR, there is a scheme-dependence in the 
extrapolation. That is, the result of an extrapolation will be dependent on the choice 
of regulator. To demonstrate this, consider the infinite- volume extrapolation of the CP- 
PACS data as the regularization scale is changed. Figure 1 shows that the extrapolations 
diverge in the chiral regime, where the non-analytic behaviour of the loop integrals 
dominate. 



Identification of the intrinsic scale: a toy model 

Is there a 'best' regularization scale indicated by lattice QCD data? By constructing 
a toy model, a method is developed to identify an in-built scale by analyzing the 
low energy coefficients. Ideal pseudodata is created with a known scale A^, and an 
examination of the renormalization flow of the low energy coefficients c; allows one 
to recover the value of the scale A^. 
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FIGURE 1. Left: example dipole extrapolation based on CP-PACS data [9], box size: 2.3 -2.8 fm. 
Right: close zoom of the regulator-dependence of the extrapolation. 



0.92 



0.91 



0.90 



0.89 



1 1 1 1 


1 1 1 1 1 1 

™-„m.' = ^^^^ 




™rZ? = <^eV= 


"\ 


m™^= = 0.5 GeV= 


\ 




~ ^ \ 
























1 1 1 1 






3.3 



3.2 



3.1 



3.0 



2.9 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 
A (GeV) 



2.8 



III 


1 1 1 1 1 1 
/ ^ ^ 


- / 
/ '1 

I \ 1 M 1 


/:i — : 

m"'""^ = 0.25 GeV= " 

JT.raax 

= 





0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 
A (GeV) 



FIGURE 2. Left: Behaviour of cq vs. A. Right: Behaviour of C2 vs. A. Both plots are based on infinite- 
volume pseudodata created with a dipole regulator at Ac = 1 .0 GeV (based on lightest four data points 
from PACS-CS). Each curve follows from fits of Eq. (1) to pseudodata, with a different upper value of 



pion mass considered. 



Based on PACS-CS data, a dense and precise set of pseudodata is generated at infinite 
volume and = 1.0 GeV, entirely within the PCR. As the pseudodata are extended 
beyond the PCR (by increasing max)' ^ scale-dependence in the renormalization is 
realized, as shown for cq and C2 in Figure 2. 

However, the correct values of cq and C2 (known by construction) are recovered at 
exactly A = Ac, where the renormalization is independent of the range of m| used. It 
is not correct in general to assume that the asymptotic value of c, at large A is the best 
estimate of the true c,-. Note that the intersection point is at the same value of A for 
both Co and C2. This indicates that the optimal regulator for the pseudodata is the value 
A = Ac. 
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FIGURE 3. Left: Behaviour of cq vs. A, based on CP-PACS data, to order ^(m^). Right: Behaviour of 
Cq vs. a, based on QCDSF data, to order ff{mji). In both plots a dipole regulator is used. A few points are 
selected to indicate the general size of the statistical error bars. 



Lattice simulation results 

In the toy model, a method was developed for extracting the intrinsic scale from some 
ideal pseudodata. This leads one to ponder if actual lattice simulation results contain an 
intrinsic scale. Since lattice data usable for an extrapolation generally extend outside the 
PGR, an optimal regularization scale may be recovered. Renormalization flow curves 
can be constructed for lattice simulation results by constraining the fit window with 
varying values of m| j^^^, sequentially including data points extending further outside 
the PGR. 

The results for the renormalization of the low energy coefficients cq, the nucleon mass 
in the chiral limit, and Cq, the magnetic moment of the nucleon in the chiral limit, are 
shown in Figure 3. Lattice simulation results from GP-PAGS are used to analyze the 
nucleon mass, and preliminary results from QGDSF are used to analyze the magnetic 
moment. Similar to the pseudodata case, a greater degree of regulator-dependence is 
observed as more data are included in the fit. There is also a reasonably well defined 
regularization scale at which the renormalization of cq is insensitive to the truncation of 
the data. This indicates that there exists an optimal regularization scale embedded in the 
lattice simulation results themselves. 

It is also intriguing that the optimal regularization scale obtained from independent 
lattice results for different observables is approximately the same value for the same 
regulator functional form. This suggests that the optimal regularization scale indicates 
the existence of an intrinsic scale in the nucleon-pion interaction. 

A measure of the systematic uncertainty in the optimal regularization scale, evident in 
Figure 3, can be obtained by constructing a ;f;^-style analysis. By plotting against the 
regularization scale A, a measure of the spread of the renormalization flow curves can 
be calculated, and the intersection point obtained. The corresponding (per degree of 
freedom) plots are shown in Figure 4. 
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FIGURE 4. Left: Behaviour of xjof 

vs. A, for the renormalization flow of cq to order 0{m]^), based on 
CP-PACS data. Right: Behaviour of Xdof renormahzation flow of Cq to order £^{m\), based 

on QCDSF data. In both plots a dipole regulator is used. 



SUMMARY 

Using finite-range regularized chiral effective field theory, a technique for obtaining an 
optimal regularization scale from lattice results was demonstrated for pseudodata, in 
addition to actual simulation results for the nucleon mass and simulation results for the 
nucleon magnetic moment. By analyzing the renormalization flow of the low energy 
coefficients with respect to A, a step-wise extension of the data beyond the PCR allowed 
the identification of an optimal regularization scale. This scale of approximately 1 .2 GeV 
is the value at which the renormalization is least sensitive to truncation of the lattice data. 
Lattice simulation results for the nucleon mass and nucleon magnetic moment yielded 
consistent values for the optimal scale. Thus an intrinsic scale was uncovered, which 
characterizes the energy scale of the pion dressing of the nucleon. 
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